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ABSTRACT 

We describe our methods for measuring the Hubble constant from Ryle Tele- 
scope (RT) interferometric observations of the Sunyaev-Zel'dovich (SZ) effect from 
a galaxy cluster and observation of the cluster X-ray emission. We analyse the error 
budget in this method: as well as radio and X-ray random errors, we consider the 
effects of clumping and temperature differences in the cluster gas, of the kinetic SZ 
effect, of bremsstrahlung emission at radio wavelengths, of the gravitational lensing 
of background radio sources and of primary calibration error. Using RT, ASCA and 
i?OSAT observations of the Abell 1413, we find that random errors dominate over 
systematic ones, and estimate Hq — 57^ kms^'Mpc^' for a an D.m — 1.0, Ha = 0.0 
cosmology. 
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1 INTRODUCTION 

Absolute distance measurements in astronomy are extremely valuable and are very difficult to 
make. The extragalactic distance scale in particular has profound implications for cosmology, 
but the distance-ladder methods have generated long-standing controversies. Direct methods of 
measuring cosmological distances, for example via gravitational lensing effects, offer alternatives 
to the distance ladder. One such method exists through combining X-ray observations of clusters 
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of galaxies with the Sunyaev-Zel'dovich (SZ) effect ( JSunyaev & Zel'dovich 197^ . This method 
takes advantage of the fact that the SZ effect depends on a different combination of the intrinsic 
properties of the cluster gas (density, temperature, physical size) than the observed X-ray flux: 
hence the physical size of the cluster, and via the angular size, its distance, can be measured. To 
see how the distance depends on the measurable quantities, consider a blob of gas of size /, electron 
number density He, at temperature Tg. The observed X-ray flux density is 

Fx'^nlf{T,)l^Df 

where Di is the luminosity distance and f{Te) is the temperature dependence of the X-ray emis- 
sivity (including Gaunt factor). The SZ effect from the same gas is 

ATsz °^ rieTel. 

If the blob subtends an angle then 

= I /Da, 

where the angular size distance = D/(l +z)^^. Rearranging, 

TT 1 ^l^X 



Ar2/(7;)0' 

where we have dropped the explicit z dependence. Hq can therefore be determined for a cluster of 
known redshift provided we can measure the gas temperature, the X-ray flux and the SZ effect, 
and provided we can model the distribution of gas in the cluster. 

This method (e.g. Cavaliere et al ( |1979P ) was suggested shortly after the SZ itself was pro- 
posed, and since the first successful detections of the SZ effect in the mid-1980s several au- 
thors have published determinations of Hq based on it, using both single-dish (e.g. (Birkinshaw 
& Hughes 1994; [Myers et al. 1997P ) and interferometric (e.g. ( Pones 1995| ; [Reese et al. 2UUq )) SZ 
measurements. Despite the advantages of the SZ method in by-passing the entire distance ladder, 
it does have its own sources of error, which must be carefully addressed. These include: 

(i) Beam switching that does not go fully outside the cluster; contamination by radio sources in 
the cluster or (for single-dish SZ measurements) in the reference beams; resolving out most of the 
total flux of the cluster (for interferometric measurements); poor spatial sampling of the cluster. 

(ii) Errors in the temperature measurement of the cluster, due to poor photon statistics, temper- 
ature substructure and clumping of the cluster gas 

(iii) Errors due to fitting a simplified smooth model of the gas distribution to a real complicated 
cluster 

(iv) Biases due to selection effects in the chosen cluster sample. 
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Measuring the Hubble Constant 3 

In this paper we describe our analysis techniques for determining Hq by combining SZ obser- 
vations from the Ryle Telescope (RT) with X-ray data from ROSAT and ASCA. We apply this 
to d ata for the cluster A1413, for which we have previously published SZ results ( Grainge et al. 
1996). Finally, we consider the various sources of error, both random and systematic, which affect 
our estimate. 



2 CLUSTER MODELLING 

We use our data to measure the Hubble constant by creating a numerical model of the temperature 
and density distributions in a cluster deduced from X-ray data, and then using these to generate 
mock RT data that can then be compared with real observations. The gas density is modelled as 
a King profile in one or three dimensions (depending on whether the cluster shows significant 
asphericity), and arbitrary temperature profiles can be specified. The cosmological parameters Hq, 
Q.M and ^2a, upon which the relative scaling of the X-ray and SZ data depend, can be adjusted 
until both sets of mock data agree with the observations. The method is essentially similar to those 
used in other SZ measurements of Hq, especially other interferometric measurements (R eese et al. 

Unlike Reese et al we do not use the SZ data to constrain the shape parameters, as the X-ray 
image contains far more information and would dominate any joint fit; we do however allow for a 
non- spherical cluster. 

2.1 Physical model 



Our modelling is based on the analytical King approximation dCavaliere & Fusco-Femiano 1976| ) 
for the potential due to a self-gravitating isothermal sphere, which gives results consistent with 
X-ray images, at least in the inner parts of the gas distribution. The cluster total density p at radius 
r from the cluster centre is given by 



p(r)=po(l + (r/r,)^ 



where rc-is the core radius of the cluster. Assuming that the gas density, p^ is related to the total 
mass density by Pg °^ pP, then 

p,(r)=p,o(l + (r/r,)') . (1) 

The average value of (3 determined by fits to X-ray surface brightness of a large number of clusters 
is found to be < |3/,y >= 0.67 (|Jones & Forman 1984j ), but we fit for P separately for each cluster. 



Although we typically use isothermal models for fitting X-ray and SZ data, we can vary the 
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temperature to investigate specific effects. We do not vary tlie X-ray emissivity constant in tliese 
models; liowever, for typical rich cluster temperatures (7-10 keV) this constant is a weak function 
of temperature when observing in the ROSAT 0.5-2 keV band. 

2.2 Simulating observations 

We generate two-dimensional images, of typically 256 x 256 pixels, of the model skies in SZ and 
in X-rays. This results in a spatial resolution in the model of ~ 50 kpc and a total simulation size 
of ~ 10 Mpc - numerical experiments show the results to be insensitive to these parameters. We 
can use either a spherical King model (equation |l]) for the density, or a 3-d modification of it in 
which core radii (r^x, /"cy, ^cz) are specified in three perpendicular directions, the directions relative 
to the plane of the sky being specified by three Euler angles. The density at a position {x,y,z) in 
these coordinates is thus given by 

We integrate the functions shown below along the rows of a three-dimensional array to generate 
the SZ brightness temperature and X-ray images: 

Ar(/,j) j,fc)r(/,j,fc)^^^L, (2) 

4 V 

X e ( , X 



hv 



X{i,j)=Y,Kn^{iJ,k)te.pL\ (3) 

k 

where Tq is the CMB temperature, tgxp is the exposure time of the X-ray image and ^ is a tele- 
scope, redshift and cluster-temperature dependent constant determining the received count rate, 
incorporating the bremsstrahlung and line emissivities, hydrogen absorption, detector bandpass 
and K-correction, expressed as ROSAT PSPC counts per second due to a volume of Im^ of gas, of 
electron density Im^-^, at a luminosity distance of 1 Mpc. L is the size in metres corresponding to 
the angular extent of each pixel Gee//; this is calculated using the cluster redshift, z, and assumed 
values of Hq, Q.m and 

We have checked that we obtain results consistent with the analytic expressions for both the 
X-ray and SZ profiles for the spherical isothermal King model. We choose not to use the analytic 
expressions in our code since the approach described above can more easily be modified to allow 
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different distributions of gas density and temperature. For example, we have been able to investi- 
gate the effect of non-isothermal temperature profiles, and also central cooling flow regions. The 
latter have no effect on the SZ signal, however, since cooling flows are considered to be in pressure 
equilibrium with the rest of the cluster. 

The brightness temperature decrement map is converted to flux density, multiplied by the RT 
primary beam response, and Fast Fourier Transformed to obtain the simulated aperture-plane re- 
sponse of the RT to the SZ decrement. The best-fit cluster parameters found for the cluster A1413 



(see section |3J|) give a profile in the aperture plane for different observing baselines shown in Fig- 
ure H]. Mock data that correspond to the real RT observations can then be generated by sampling 
this array at the observed coordinates in the aperture plane. 

We also simulate the process of observing the sky with the ROSAT PSPC or HRI. The model 
X-ray sky is convolved with the X-ray telescope point-spread function, and a suitable background 
level added. For making model fits a smooth background is used; we can also make images in 
which each pixel value (signal plus background) is replaced with a sample from a Poisson distri- 
bution of the same mean, thus generating realistic X-ray images. 

We tested the program to ensure that it gave the expected dependence upon central number 
density, (3, core radius, gas temperature and value for Hq and to check that the results for the 
simulations were not dependent upon pixel size or overall array size. 



2.3 Predictions of SZ observations 

We have predicted the effect that changes in various cluster parameters would have on our RT 
observations. In particular we were concerned about the behaviour of the flux density measured on 
the shortest RT baselines, since these are the most sensitive to the SZ effect. (The shortest physical 
RT baseline is 870 X, but at for lower declinations this can be projected down to a minimum of 
the antenna diameter of 640 X.) The results here were obtained using the model parameters for 
Abell 1413, but are typical of all the rich clusters we observe. 

We simulated a cluster with using the best-fit parameters that we have found for A1413 (see 
section |TT|) but varied the core radius. We found that the central decrement is directly proportional 
to the core radius, as expected, and that the flux density measured on the 870 X baseline quickly 
began to flatten off (see Figure as the RT starts to resolve out the SZ signal. This resolving 
out starts to become significant for core radii of greater than about 70 arcseconds for the 870 X 
baselines. 
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In a further test we varied the central number density, hq, for various core radii to give the 
same observed flux density using a 870?i baseline, and recorded the central temperature decrement 
that these models would give. The results are also shown in Figure 0. The central decrement does 
not vary greatly with core radius, especially around 60 arcseconds which is the expected size for a 
cluster with core radius of about 200 kpc at moderate redshift. 

We have also assessed the effects of redshift and cosmology on the SZ signal seen. We sim- 
ulated observations of clusters with the same intrinsic parameters lying at different redshifts. 
First we calculated the clusters' angular core radii, corresponding to a constant physical radius 
of 250 kpc, assuming Einstein-de Sitter (fl = 1.0) and Milne {'D. = 0.0) cosmologies. The results 
for observations with baselines of 870 and 1740 X are shown in Figure ^. It can be seen that the 
RT has similar sensitivity to a cluster at redshifts between 0.2 and 10 on the shortest baseline in 
both cosmologies. On the 1740 X baseline the signal peaks for clusters at redshifts between 0.3 and 
4 in an Einstein-de Sitter cosmology, but tends to a maximum value in a Milne universe beyond 
a redshift of 1 . We therefore conclude that observations of the SZ effect with the RT should be 
possible for clusters at any redshift from z = 0. 1 to z = 10. 

Although the SZ effect predicted to be detected by the RT from a cluster of given physical size 
is largely independent of cosmology, the value of Hq that we will calculate from observations will 
be affected by the values of Q.m and Q.^_ that we adopt. The true value of Hq and the calculated 
value Hq^'^ will be related by 

r/calc ncalc 
^ D 

where Uf'^'^ and Da are the calculated and true angular distances to the cluster respectively. The 
correction factor jH^ is shown in Figure ^. 

2.4 Fitting to X-ray images 

Since KOSP^ X-ray images contain far more information about the gas distribution in a cluster 
than RT SZ data, we use the X-ray image to model the density distribution. We assume the cluster 
is a triaxial ellipsoid with two principal axes in the plane of the sky, and with the core radius along 
the line of sight equal to the geometric mean of the other two core radii. We then use a downhill 



simplex fitting method, implemented by the routine AMOEBA ( press et al. 1986D , to optimise the 



central density, core radii, ellipsoid centre, position angle and |3 parameter. The likelihoods of the 
data given the model are calculated as the product of the Poisson likelihoods of obtaining the 
observed count at each pixel in the X-ray image, given the parameters. It is necessary to use the 
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Poisson statistic rather than in order to avoid biases in the regions of low count rate where the 
two statistics differ significantly. This method also allows arbitrary regions to be excluded from 
the fit; this is used for example to exclude point sources in the X-ray image. 

We use vignetting-corrected X-ray images (i.e. divided by the image of effective exposure 
time) for the fitting. Although this results in a distortion of the photon statistics across the map 
(partly due to the presence of the non-vignetted instrumental background), this effect is small in 
the central part of the image where all the cluster signal lies. We also use images that are binned 
so as to under-sample the PSF (pixel size equal to FWHM of PSF) so that the position dependence 
of the PSF is negligible over the region of interest. 

There is a degeneracy between |3 and 0c, and n^, in that almost equally good fits can be obtained 
by changing the value of [3, and re-fitting for 0^ and no- However, although the resulting value of 
the central temperature decrement is changed for these different fits, the flux density observed 
by the interferometer is almost unchanged. This is because the observed flux density depends on 
both the central decrement and the shape of the cluster, and on our most sensitive baselines the 
two effects almost cancel. The effect of this can be seen in Figure ^ which shows contours of 
likelihood of the X-ray fit in the (3-0c plane overlaid with contours of the predicted flux density 
observed by the RT. The flux density contours (which are used to calculate Hq) lie parallel to the 
P-0C degeneracy, which thus has little effect on the value of Hq found. Reese et al. ( [2000[ ) find a 
similar result for the OVRO/BIMA arrays. 

2.5 Fitting to SZ data 

Interferometers measure the Fourier transform of the sky brightness distribution at discrete points 
defined by the orientation of the baselines, with independent noise values at each point. Invert- 
ing the data to form an image results in both signal and noise being convolved with the Fourier 
transform of the sampling function (the synthesised beam), resulting in long-range correlations 
across the map. With well-filled apertures and high signal-to-noise, this is not a problem. Since 
the RT has few baselines on which the SZ effect is detectable, these correlations are significant, 
ie the synthesised beam has large far-out sidelobes. It is therefore very difficult to do meaningful 
fitting to the data in the image plane. Accordingly, we fit our data to models in the aperture plane, 
although we also make maps, particularly in order to identify positive sources in the field. 

Once a fit has been made to the X-ray image using an assumed value of Hq, an SZ sky model 
and simulated RT aperture are generated as described above. Variations in Hq result in a simple 
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scaling of the amplitude of the SZ signal. The calibrated, source-subtracted RT visibility data are 
read in and a statistic calculated between the real and mock data as a function of amplitude 
scaling; we multiply by a prior that is uniform in log space (since Hq is a scale parameter) to find 
the peak and extent of the posterior probability distribution. 



3 CALCULATING Hq FROM A1413 

A1413, at a redshift z = 0.143, is a massive and highly luminous cluster with a ROSAT-band 
luminosity of 2.0 x 10^^ W ( [Allen et al. 1995[ ). We have previously published SZ data on the 
cluster Abell 1413 ( [Grainge et al. 1996| ) but did not derive an Hq value there due to the significant 
ellipticity of the cluster and the presence of a moderate cooling flow, both of which complicate the 
process. Here we calculate Hq from A1413 taking both these effects into consideration. As well as 
our RT SZ data, we use an observation made with the ROSAT PSPC, and a temperature derived 
from ROSAT HRI and ASCA data. 



3.1 X-ray image fitting 

We fitted an ellipsoidal model (as described in section [2^ to the ROSAT PSPC image observed 
on 1991 November 27, which has an effective exposure of 7696 s. We use only the hard (0.5- 
2 keV) data in order to minimise the effect of Galactic absorption; the image is shown in Figure 
along with the RT SZ image. Allen et al. ( [1995[ ) find that A 14 13 has a cooling flow with a mass 
deposition rate of 200 Moyr^ within a radius of 200 kpc of the centre, and that the temperature in 
this region is depressed by some 20% from the value in the outer regions. Therefore we excluded 
the X-ray data within the cooling flow radius from our analysis, along with several X-ray point 
sources which appear in the image. Keeping the cooling flow region in the fit would result in a 
more compact model (higher [3, smaller core radius) that would not reflect the true distribution of 
the gas, most of which is outside the cooling flow. 

Allen and Fabian ( |1998D use ASCA and ROSAT HRI data to determine the temperature of 
A 14 13 in two models; one assuming the gas is isothermal, and one also fitting for a cool com- 
ponent, with the fraction of the total emission from the cool gas normalised by the fitted mass 
deposition rate of the cooling flow. These methods yield 7.51q 3 keV and S.S^Qg keV respectively 
(90% confidence errors). The latter value is a better estimate of the temperature of the bulk of 
the cluster gas, outside the cooling flow region. Since we have excluded the cooling flow region 
from our fit to the surface brightness image, we use this higher temperature as the more consistent 
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model. The value of Hq that would be obtained from using the cooler temperature can be found 



from the scaling given in section |47r2| . The X-ray emissivity constant was calculated assuming an 
absorbing H column of 1.62 x 10^^ ni^^, a metallicity of 0.3 solar and a temperature of 8.5 keV, 
giving a value of 3.41 ±0.24 x 10^^^ counts s^^ from 1 m-^ of gas of electron density 1 m^-^ at a 
luminosity distance of 1 Mpc. 

The best-fitting model parameters were (3 = 0.57, core radii of 53" and 35" with a position 
angle of the major axis of 0.1° and a central density no = 1.51 x 10^/z^/^m^^ (assuming a core 
radius in the line of sight of 43 = (53 x 35)^/^ arcsec), where Hq = 100/zkm^'Mpc^^. These 
parameters agree well with the fit by Allen et al ( |1995D . The ROSAT image, our X-ray model 
and a map of the residuals are shown in Figure 0. The residual map has zero mean in the regions 
included in the fitting, although the excess flux in the cooling-flow region and a point source stand 
out. The Poisson noise also rises towards the centre of the cluster as the square root of the signal 
level. The reduced value of for this fit is 1.02; the mean log likelihood over the region of 
the fit is L = —1.585, and the expected distribution of L, as discussed in Grainge et al. ( |2001D , 
is L = — 1.585±0.021. We therefore conclude that the fit is a good one. There is a degeneracy 



between P and the core radii (see section which lies parallel to lines of constant predicted flux 



density observed by the RT. Figure ^ shows that acceptable fits to the X-ray map correspond to 
uncertainties in the predicted SZ flux density of ±14^Jy and ±24;uJy at 67% and 95% confidence. 
Therefore we estimate that the error in Hq that arises solely from our X-ray fitting procedure is 

14 
650 



2x^ = 4% 



3.2 Relativistic Corrections 

The formula for the SZ effect quoted above (eqn ^ uses non-relativistic approximations, which 
become more significant at higher gas temperatures and higher observing frequency. Challinor and 
Lasenby ( |1998D show that a better approximation in the Rayleigh- Jeans region (correct to second 
order in Te) is 

Ar ^ / nkBTe , 123 fhsT,^ 



Tq \ lOnigC^ 40 \mgC^ 

For A1413 this correction amounts to some 2.5%, and we use the corrected value in the results 
below. 
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3.3 Fitting for Hq 

The SZ data used were those described by Grainge et al ( |1996D , who also describe the procedure 
for subtracting the effects of radio sources within the field. We used only the visibilities from 
the RT baselines shorter than 2 kX for the fitting, since the SZ effect is completely resolved out 
on longer baselines. The mock data based on our X-ray model were compared with the source- 
subtracted data from our observations to find best-fit values of Hq and the corresponding value 
of uq; these were Hq = SV^jj kms^^Mpc^^ and hq = 9.0 x 10-^ m^-^ (for an Einstein-de Sitter 



universe and apply the correction for primary flux calibration discussed in section H-.lD . A plot of 
the posterior probability distribution is shown in Figure The central SZ temperature decrement 
in this model is 825 ^uK. The quoted (l-o) error on Hq is just that due to the noise on the SZ data 
and the fitting to the X-ray image; other sources of error are discussed below. 

4 OTHER ERRORS IN THE CALCULATION OF Hq 

There are other sources of error which may affect our estimate of Hq. These arise in the radio 
observations, in the X-ray data, and in the modelling process. 

4.1 Primary flux calibration of the RT 

The RT is primary flux-calibrated every day on either 3C 48 or 3C 286, assumed to have flux 
densities in the observed Stokes' parameter I+Q of 1.70 and 3.50 Jy respectively. These fluxes are 
derived from Baars et al. ( |1977[ ). We estimate that the daily random error in primary calibration is 
5%. Since the total integration time on A1413 is 65 days, the overall error in primary calibration 
due to variation during each day is ±0.6%. Since the calculated value of Hq is proportional to the 
square of the measured SZ flux, this leads to an error of ±1.2% in Hq. 

Measurements by the VLA in 'D' configuration have shown that the flux densities of both 
3C 48 and 3C 286 show slight variations with time (|VLA calibration manual 1996D. This time 



variability is small, less than 5%. Measurements in 1995, roughly contemporary with our observa- 
tions, showed that the true flux density of 3C48 was 3.4% lower than the Baars et al value, while 
that of 3C286 was 1.2% higher. Since roughly equal numbers of the individual RT observations 
of A1413 were calibrated from each source, we average the errors and take the likely systematic 
error as 1 . 1 ± 2.5%. This results in a reduction in Hq of 2.2%, and a further error term of ±5%. 
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4.2 Source subtraction 

We removed the effect of contaminating radio sources using higher resolution data from both the 
longer baselines of the RT and from VLA imaging of the field. However, the accuracy to which 
the sources can be removed depends on the level of uncalibrated amplitude and phase errors in the 
RT. From maps of secondary calibrator sources we have found that the dynamic range of the RT 
is typically 100:1. Since the total flux density of the six sources removed from the A1413 data is 
3.1 mJy, the maximum spurious flux remaining after source subtraction will be 30 /jJy. However, 
taking into account the distribution of the sources on the sky and the contribution each makes at the 
position of the SZ decrement, we find that the maximum unsubtracted contribution of the detected 
sources is only 3 ^uJy. 

However, a more significant contribution is from the sources which were not subtracted due 
to their being too faint to detect. To calculate the confusion noise contribution from unsubtracted 
sources, we take the 8.4 GHz /jjy source counts of Windhorst et al ( |1993[ ), and extrapolate to 
15 GHz using an effective spectral index calculated from their 5-8.4 GHz spectral measurements. 
Integrating the source counts over the beam area corresponding to our shortest baseline, from zero 
flux density up to our source detection limit of 120 /jJy, and taking into account the increased 
source population in clusters compared to the field ( [Cooray 1999[ ), we find a residual confusion 
noise of 60 /jJy, which adds in quadrature with the thermal noise. However, since we estimate our 
thermal noise from the scatter in the visibilities, the confusion noise is already included in our 
noise estimate, and does not need to be added in again. 

4.3 X-ray emission constant 

The X-ray emission constant, K, used in simulating the X-ray data is dependent upon the telescope 
detector response in the appropriate frequency band, the K-correction for redshift, the cluster tem- 
perature and the absorption due to Galactic hydrogen. We estimate that this constant could be in 
error by ±7% leading to a ±7% error in Hq. 

4.4 Extended radio emission 

Some clusters of galaxies are known to have diffuse ('halo') emission on arcminute scales. Any 
such emission at 15 GHz would be resolved out on the longer RT baselines used for source sub- 
traction, and would therefore contaminate the SZ data. Observations at 1 .4 GHz from the NVSS 
and FIRST surveys show that there is 1.9 mJy of flux associated with the central cluster galaxy 
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which is extended on scales of ~ 45", but none on larger scales. This source was detected (and 
indeed slightly resolved) in the RT longer baselines and was removed from the data (source 3 in 
Grainge et al ( |1996D ). The spectral index of this source u}^^ is 0.8, typical for a cluster galaxy 
but quite atypical for a diffuse halo source, where one would expect > 1.4 ( Panisch 1982| ). 
Also, halo emission is unknown in clusters with strong cooling flows such as A1413. We therefore 
conclude that there is no evidence for any unsubtracted diffuse emission in A1413. 

4.5 X-ray estimate of gas temperature 

The temperature we use has a quoted 90% confidence error of 1^^%. Converting this to a l-o 
errors for consistency gives which will lead to a error in Hq. 

4.6 Cluster ellipticity 

X-ray images of clusters show that clusters are typically elliptical with ellipticities of up to 1.2:1 
being common. Therefore one cluster is not sufficient to provide a definitive value of Hq. Making 
the necessary corrections to equation [l|, one finds that the calculated value //q^'*^ is related to the 
true value by: 

r/calc ^ u true^-L 
7|i"' 

for a cluster with line-of-sight depth and observed diameter /|| . Therefore the derived Hq from a 
sample where there is no bias in the orientation of clusters will form a distribution with a geometric 
mean close to the true value. It is thus necessary to observe such a sample of clusters in order to 
reduce the uncertainty due to this effect. An X-ray surface-brightness limited sample will be biased 
towards selecting clusters elongated towards us. However, an X-ray sample selected on X-ray flux 
will be independent of cluster orientation. 

In a case such as the present where one is estimating Hq from a single cluster, it is appropriate 
to make some allowance for orientation uncertainty. For A1413 we have assumed that the line- 
of-sight depth through the cluster is the mean of the two elliptical axes in the plane of the sky. 
This is an unbiased estimator of the true depth, assuming the cluster is drawn from an unbiased 
population ( [Grainger et al. 200Tt ). (Sulkanen ( |1999D finds a similar result taking the arithmetic 



mean.) In order to estimate the error introduced by this assumption we compare model distributions 
of elliptical clusters with those observed in a sample of ROSAT-selected clusters over a redshift 
range of z = 0.1-0.5. We find that the effect of cluster ellipticity adds an error of 14% for each 
cluster to the calculated value of Hq. 
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4.7 Effect of temperature structure 

Considerat ions of hydrostatic equilibrium, as well as observations on nearby clusters ( Markevitch 
et al. 1998) show that the isothermal assumption breaks down at radii much bigger than Vc, the 
temperature falling by a factor of about 2 at 6rc. This can in principle have a significant effect on 
the derived value of Hq, tending to bias it downwards ( p^nagaki et al. 1995| ; [Roettiger et al. 1997| ). 



For A1413, one core diameter corresponds to 2 arcminutes, which is approximately the scale to 
which the RT is sensitive. Gas on larger scales than this does not affect the flux density seen by 
the RT, despite the effect it has on the SZ central decrement. ASCA observations do encompass 
larger scales, but the temperature measurement is necessarily emission-weighted, and hence is a 
good measurement of the temperature of the gas that the RT sees. We therefore conclude that SZ 
observations by the RT are not going to be significantly affected by the presence of a cold halo of 
gas, and so we can calculate Hq without requiring knowledge of the gas temperature at large radii. 
A future paper will deal with this effect in more detail. 

4.8 Clumping of the intracluster gas 

X-ray imaging on nearby, large angular size clusters has found no evidence for clumping, and has 
placed constraints on the degree of clumping which may be present ( pabian et al. 1994] ). However 
it has been suggested that the scatter in the temperature-luminosity correlation for clusters is due 
to the existence of clumping below this level and further that the degree of clumping is correlated 



with the strength of the cooling flow ( |Kim et al. 1991D . If the clumps are close to being in pressure 



equilibrium with their surroundings the clumping will have no effect on the SZ signal. The X-ray 
luminosity will increase but the effect of this in the Hq estimate will largely be cancelled due to 
the decreased emission-weighted temperature. The net result would be a modest underestimate of 
Hq. Initial simulations indicate that fractionation leading to an underestimate of 2% in Hq would 
be consistent with the temperatures calculated by Allen and Fabian ( |1998D . 



4.9 The kinetic SZ effect 

Bulk motion of the cluster gas along the line of sight will result in a further distortion of the CMB 
spectrum due to the Doppler effect. Watkins ( |1997| ) finds that the 1-d rms peculiar velocity of 



clusters is Oy, = 265^^^ km s^^ (at 90% confidence). At 15 GHz, a cluster with a peculiar velocity 
of 265 km s^ will have a kinetic SZ effect with a magnitude 2.5% of that due to the thermal 
SZ effect and so introduce a 5% error in the determination of Hq from a single cluster. Since 
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randomly selected clusters will have random directions of peculiar motion, the error introduced 
into Hq calculations can be reduced by averaging over a large sample of clusters. 



4.10 Radio bremsstrahlung 

It is possible that the electron cluster gas could emit sufficient thermal bremsstrahlung radiation at 
15 GHz to significantly affect the SZ decrement which we observe. We calculated the emissivity 
of the hot thermal component of the gas using our cluster model, taking the Gaunt factor to be 
10 at a frequency of 15 GHz. We find that the total integrated flux density from the entire cluster 
is 15 ;uJy. This emission will be distributed like the X-ray emission, and the resolving-out of our 
shortest baseline means that we will observe only about 20% of it. Since the emissivity due to 
bremsstrahlung is proportional to rp-, it is possible that the effect of thermal bremsstrahlung will 
be much greater in a cluster which has a strong cooling flow and thus a high gas density at the 
cluster centre ( P^arter 1978| ; [Schlickaiser 1991D . We investigated this possibility by using a good 



empirical fit to the form of a cooling flow, determined from X-ray observations, allowing the 
density to increase as 1 /r for r within the cooling flow radius of 100-200 kpc until r = 10 kpc. We 
find that this cooling flow contributes a further 4 /jjy of flux density. Therefore we conclude that 
even for this worst case example, and assuming that none of the cooling-flow flux is resolved out 
by the RT, the effect of thermal bremsstrahlung will be to overestimate Hq by up to 2%. 

4.11 Effects of gravitational lensing 

The strong gravitational potential of the cluster will inevitably lead to lensing of background radio 
sources in the field. Refregier and Loeb ( |1997| ) have calculated the systematic bias this introduces 



into Hq calculations. Assuming P = 0.6, using a 2' beam and subtracting sources above 100 //Jy, 
we find that approximately 0.5% of the SZ temperature decrement may be due to this effect. This 
leads to an underestimate in our value to //q of 1%. 

4.12 Combining the errors 

All the errors listed above are summarised in table [l|. Writing the result as a function of the sys- 
tematic quantities concerned, we find 

(1.64 Jy , 3.54 Jy\ 2 ,^ _ fiox 

xf:^Vxf4^xri +0.085 



8.5keV/ V57"/ V ' 1000 km s-i 
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(1 "^brems \ f ^ , "^lens \ , — In* —1 
1 X IH km s Mpc 
640A/Jyy V 640A/Jyy 

where 53C48 and 53C286 are the true primary calibrator flux densities (I+Q), K is the X-ray emissiv- 
ity constant as defined in section |TT], Te is the gas temperature, 0^ the line-of-sight core radius, 
the line-of-sight peculiar velocity, and Sbrems and 5iens are the flux densities due to bremsstrahlung 
emission and lensed background sources respectively. Inserting our best estimates for the errors in 
these quantities, we find in the worst case, where the errors conspire to add together linearly, 
Hq = 57^27 kms^^Mpc^^. However, since these errors are all independent of each other, we 
can justifiably combine them in quadrature with the random error. Doing this we obtain Hq = 
57^j^ kms^^Mpc^^ The largest errors are due to the uncertainty in X-ray temperature, the noise 
on the SZ measurement and uncertainty in cluster line-of-sight depth. 

5 CONCLUSIONS 

We have described our methods of determining the Hubble constant from SZ and X-ray observa- 
tions of clusters of galaxies, giving careful consideration to sources of both random and systematic 
error. We have used observations of the cluster A1413 to estimate Hq. We conclude that: 

(i) The error due to unsubtracted radiosources is comparable to the thermal noise, and adds in 
quadrature with it; our estimate of the noise from the scatter in the visibilities includes both effects. 

(ii) Fitting the X-ray emission using a King model gives a degeneracy in the |3, plane. The 
predicted flux density that the RT will observe based on any of these good-fit models is almost 
constant, with the result that the error in Hq that is introduced solely by our X-ray fitting is small, 
approximately ±4%. 

(iii) There is a random uncertainty in the calculation of the X-ray emissivity of the cluster gas 
which we estimate to be 7%. 

(iv) The orientation of any non- spherical cluster causes a random uncertainty in the value of Hq 
of some 14%. This error may be reduced by observing a complete sample of randomly-oriented 
clusters. 

(v) A colder cluster atmosphere at large radius will not affect our determination of Hq much, 
but clumps of cool gas within the cluster may result in an underestimate the value of Hq of up to 
2%. 

(vi) The kinetic SZ effect will cause a random uncertainty of approximately 2.5% the value of 



© 0000 RAS, MNRAS 000, 000-000 



16 K. Grainge et al 

the thermal SZ effect. This error may be reduced by observing a complete sample of randomly- 
oriented clusters. 

(vii) We estimate that the contribution of thermal bremsstrahlung at radio wavelengths will 
cause an overestimate of Hq of 2% in the worst case. 

(viii) Gravitational lensing will cause the over-subtraction of flux from contaminating back- 
ground radiosources; this is a small effect causing an underestimate in Hq of about 1%. 

(ix) The dominant sources of error are thus the thermal/confusion noise in the SZ measurement, 
the temperature of the gas, and the unknown line-of -sight depth. Combining these sources of 
random and systematic error in quadrature, we determine from our observations of A1413 a value 
ofHo = Sltjl kms-^Mpc-i. 
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Source of error 


% error in Hq 


Noise on SZ 


±28% 


Fitting to X-ray map 


±4% 


RT primary calibration 


±5% 


Source confusion 


±16% 


X-ray emission constant 


±7% 


Cluster gas temperature 


+20/ -12% 


Ellipiticity 


±14% 


Clumping 


+2% 


Kinetic SZ 


±5% 


Radio bremsstrahlung 


-2% 


Gravitational lensing of background sources 


+1% 



Table 1. Error budget for Ha determination from A1413. For each source of error a + sign indicates that the true value of Ho may be more than the 
calculated value. 
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Figure 1. Simulated S Z fl ux as a function of projected RT baseline. The cluster model used for this simulation was the best fit to the ROSAT map 
of A1413 (see section [ij] ). During the course of a 12 hour run the RT observes with projected baselines b > 640X; baselines shoiler than this are 
shadowed by other antennas and the data are flagged. 
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Flux Against Core Radius, Gas Density Constant Central Decrement Against Core Radius, Flux Constant 




Core Radius (arcsec) 



Figure 2. (Left) Observed flux density on a baseline of 870 X against cluster core radius, with constant central gas density and temperature. (Right) 
Central decrement against core radius, with constant observed flux density of 550 /Jy on a baseline of 870 X, showing the relative insensitivity of 
observed flux density to central decrement at this resolution. 
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Figure 4. Correction factor to calculated value of i?o for different cosmologies, relative to an Einstein-de Sitter model. 
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Figure 5. Greyscale and solid contours: likelihood of fit to X-ray data, as a function of P and core radius of the major axis, marginaUsed over no but 
with the axial ratio and position angle (which are very well determined) fixed. Contours are 67% and 95% confidence. Dotted contours: predicted 
mean flux density on the shortest RT baseline, for h = 0.5. The degeneracy between P and the core radius clearly has little effect on the predicted 
flux density, and hence on the fitted value of Hq. 
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Figure 6. X-ray and SZ images of A1413, from Grainge et al (1996). The greyscale is the ROSAT PSPC image; it runs from 2 to 80 counts in 
a 15" X 15" pixel. The contours show the CLEANed map of the naturally weighted RT 0-1 kX visibilities with restoring beam 169" x 110" at 
position angle 2.7°; contour levels at -450, -375, -300, -225, -150, -75, 75 ijjy. 
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Figure 7. PSPC image, our X-ray model, and the residuals from the fit. Contours are every 1.0 counts per 8" pixel. The cooling flow and an X-ray 
point source which were not included in the fitting procedure are evident in the residual map. 
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Figure 8. Posterior probability plot for values of Hq from the fit of SZ decrement to the X-ray derived cluster model. The width of the plot indicates 
only the errors from the visibility data. 
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